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Abstract 

We propose a method for computation of stable and unstable sets 
associated to hyperbolic equilibria of nonautonomous ODEs and for com¬ 
putation of specific type of connecting orbits in nonautonomous singular 
ODEs. We apply the method to a certain a singular nonautonomous real 
Ginzburg-Landau type equation, which that arises from the problem of 
formation of spots in the Swift-Hohenberg equation. 
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1 Introduction. 

The aim of this paper is to present a method for computation of connecting 
orbits between hyperbolic equilibria in nonautonomous ODEs. Although the 
algorithm we propose is general and may be applied to an ODE in any (finite) 
dimension, the main motivation for undertaking this research was the problem 
of finding transversal solutions for the following nonautonomous real Ginzburg- 
Landau type equation 

Arr + ~ ~ ^2 = A ~ a3 ’ (k 
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where A : [0, oo) —> K. satisfies 


T(0) = 0, lim A{r) = 0. 

r—>• oo 


( 2 ) 


This question comes from McCalla and Sandstede work IMS] on spots in the 
Swift-Hohenberg equation. The existence of solutions ([!]) satisfying ([2]) is listed 
there as an assumption in the construction of localized solutions called the B 
spots. The work of Scheel [3] (references given there) indicates that bounded 
solutions of problem ( TJ2 ) might be also interpreted as rotating wave solutions 
to reaction diffusion systems in the plane. 

It was pointed out in [MS] that the work by Scheel |S| on (l|2), where the 
existence of an infinite number of geometrically different solutions has been 
proved, contains a gap. Recently, a computer assisted proof of the existence of 
one such solution was given in |BOW| . The following theorem summarizes our 
main results about the system 0 . 


Theorem 1 For n = 1, 

such that: 


,6 there exists a solution A n : (0,oo) 


to [ 7p ] 


• lim r _ >00 (A n (r),T' n (r)) = (0,0), 

• \im. r ^ 0 +{A n {r),A' n (r)r) = (0,0), 

• The function r —> A n (r) has exactly n local extremes in the domain r £ 
(0, oo). 

The solution A n is transversal in the following sense. The sets 

W u = {(r 0 ,4), O) £(0,oc)xfxR 

the backward solution of 0 A(r) with initial condition 
A(r 0 ) = A 0 , A'{r 0 ) = A' 0 satisfies lim r _ >0 (^(?'), A!(r)r) = 0} 

IT S = {(r 0 ,A 0 ,A' 0 )) e(0,oc)xf xR 

the forward solution of 0 A(r) with initial condition 
A(r 0 ) = A 0 , A'(r 0 ) = A' 0 satisfies lim r _ ) . 00 (T(r), A!{r)) = 0} 


are two dimensional immersed manifolds in the extended phase space (r,A,A') 
and they intersect transversally along the solution curve defined by A n (r), i.e. 
(r, A n (r),A' n {r)). 


Plots of the six connecting orbits resulting from Theorem [I] are shown in 
Fig- H The solutions differ in the number of local extremes or the number of 
sign changes. The orbit for n = 1, with one maximum only, is the one already 
proved in |BGW] . We present a proof of the existence of six such orbits, but 
without much difficulty it would be possible to validate the existence of several 
orbits with more local extremes. The question of the existence of an infinite 
number of geometrically different solutions for © is still open. 
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It should be noted that Ryder m proved the existence of infinitely many 
geometrically different solutions to a similar problem 

2 

u" H— u' — u + u 3 = 0. (3) 

r 

plus boundary conditions analogous to ([2|, which arises as the equation for the 
radial solutions in the Ginzburg-Landau equation in 3D. Apparently, Ryder did 
not proved the transversality. 

The Ryder’s method applies to a wide class of equations, however due to the 
presence of term , equation 0 does not belong to this class. 

Using our approach we can obtain several first solutions also for § (ap¬ 
parently without much difficulty because it appears to be numerically easier 
than 0), but the number of solutions will be always finite in contrast with the 
Ryder’s approach. Our shortcoming is due to the following: the lack of mathe¬ 
matical understanding was replaced by the computer power. What is missing is 
the understanding how the whole (un)stable set behaves away from the origin, 
and it is exactly the term ^, which causes the biggest problem. 

Our approach is based on the shooting method (see for example [HM| ), which 
is a very well known technique in the dynamics of ODEs and boundary value 
problems for ODEs. 

To realize the shooting method for problem 
following issues 

1. the singularity for r = 0, 

2. the estimation of stable manifold of the point (A, A') = (0,0) for non- 

autonomous equation 0- 

Issue 1 is dealt with via the change of the independent variable p = lnr. Then 
the equation 0 becomes 

w ” — w/ 4 = e 2p (w — w 3 ). (4) 

Now we seek a solution such that liuip^-oo w(p) = 0, so we need estimates for 
the unstable set of (w = 0, w' = 0) for p —>• —oo for our non-autonomous ODE 
Q, which is basically issue 2 above. 

Issue 2 in our work is split in two separate problems 

• obtaining analytically bounds for the stable manifold in the neighborhood 
of the hyperbolic fixed point, which are valid up to some ’macroscopic’ 
distance from it, 

• the globalization of the local (un)stable manifolds using the rigorous nu¬ 
merical integrator from CAPD CAI'Dl . 

Obtaining explicit bounds on the local (un)stable manifolds associated to 
hyperbolic equilibria in nonautonomous systems is an interesting issue itself 
and most likely this is the main mathematical contribution of this paper. For 


ip i we have to deal with the 
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Figure 1: The six transverse connecting orbits resulting from Theorem [TJ In 
two upper rows in red we have plotted the separatrices (stable and unstable 
manifolds of (0,0)) of system A rr = A — A 3 , which is the formal limit of ([l]) for 
r —>• oo. 
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this end we adopt the method of cone conditions based on quadratic forms from 
[ZCC, iSZ^ to the case of nonautonomous systems. The main novelty of this 
paper in this context is Lemma [4] which gives explicit bounds on the derivative 
of the parameterization of the stable manifold of a hyperbolic equilibrium with 
respect to the time variable. The method is general and may be applied to 
any smooth nonautonomous system in arbitrary finite dimension. We believe it 
should be an important ingredient when attacking the question of the existence 
of infinite number of geometrically different solutions to 01 which is the next 
logical step for this project. 

To demonstrate the effectiveness of our estimates of local (un)stable mani¬ 
folds in Sections [4] and [5] we present easy calculations showing that using our 
approach we can estimate invariant manifolds for the system Q up to the dis¬ 
tance of order 0.1 from the origin. Outside of this region the stable and unstable 
manifolds of (flj) and p} are continued to obtain a connecting orbit using a rig¬ 
orous ODE integrator from the CAPD library jCAPPj . but for this to succeed 
it is essential that our local bounds are valid up to a ’macroscopic’ distance 
from the origin, otherwise the execution of the rigorous numerics will take very 
long time and the quality of bounds will be poor. Let us emphasize that to glue 
unstable and stable manifolds we need to integrate the system and associated 
variational equations over a finite time interval, only. Our local estimates on 
the invariant manifolds and bounds coming from the rigorous integration were 
sharp enough to fulfill the assumptions of the shooting method that guarantees 
the existence of six connecting orbits as stated in Theorem [T] 

Let us comment about the computer assisted proof from IBGWj . There the 
problem is formulated in the functional analytic language. Differential equation 
0 with the help of Green function is replaced by an integral equation on finite 
domain and conditions ([2]) is encoded in the choice of the function space. There¬ 
fore obtaining a desired solution of 01 is reduced to the fixed point problem 
for some nonlinear integral operator in some infinite dimensional function space. 
This should be contrasted with our approach, where we work in the phase space 
of our ODE. 

The content of this paper can be described as follows. In Section [2] we recall 
the main geometric tools used in our approach to the (un)stable manifolds. In 
Section[3]we present theorems which provide us with bounds on the derivatives of 
(un)stable manifolds. In Sections [4] and [5] we give explicit bounds for (un)stable 
manifolds and theirs derivatives for 0 and Q. In Sectionjbjwe give a computer 
assisted proof of Theorem [l] 


2 Geometric tools for invariant manifolds of fixed 
points 

To make this paper reasonably self-contained in this section we gather all nec¬ 
essary definitions and geometric theorems from |ZGi[ IWZl j ZCCl IWZ1 related 
to the rigorous investigation of invariant manifolds of fixed points, with rnodifi- 
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cations required for the setting of non-autonomous ODEs. 


2.1 Notation 

For a matrix A £ R” xfc by A T we will denote its transpose. Let Q be a square 
matrix. With some abuse of notation we will denote by the same letter a 
symmetric matrix Q and associated quadratic form Q(x ) = x T ■ Q ■ x. To 
avoid ambiguity we will always use brackets Q(-) to indicate that Q is treat 
as a quadratic form and Qx or Q ■ {x\ — X 2 ) when we compute matrix-vector 
product. For a bilinear map M : V x W —t R we will use the norm ||W|| = 
max |M|<i,|M|<i \M(v,w)\. 

For u £ N, p £ R“ and r > 0 by B u (p,r) we will denote an open ball 
of radius r centered at p. We will also often use B u = B u { 0,1). For a map 
/i:[0,l]xI->X and t £ [0,1] we set h t {x) = h(t, x). 

For a function f(t,x) we set D t f(t,x) = |£( t,x ) and D x (t,x) = S£(t,x). 

2.2 Horizontal and vertical discs 

The main tools for proving the existence of local (un)stable manifolds in a given 
domain are isolating blocks and the cone conditions. The invariant manifolds 
are then constructed as graphs of smooth functions in some isolating blocks, as 
horizontal/vertical discs in h-sets. The size of an h-set gives an explicit range 
in which we can parameterize the local invariant manifold. The cone conditions 
are used to obtain both the existence and the local uniqueness of the invariant 
manifold as well as bounds on the derivatives of this parameterization. 

Definition 1 fZUU Def. 1] 

An h -set N is a quadruple (\N\,u(N), s(N),Cn) such that 

• \N\ is a compact subset of K”, 

• u(N ), s(N) £ {0,1,2,...} are such that u(N ) + s(N) = n, 

• cat : R 71 —► R ra = x R 5 ^ is a homeomorphism such that cjv(|2V|) = 

B u {n) x ^s(AT)- 

We set 


dim (N) 

= n. _ 

N c 

= -Bu(jv) x -B s (ao, 

Nf 

= dB u ( jy) x £? s ( at) 

N+ 

= AT) x <9f? s (AT) 

N~ 

= c n\N-), 

N+ 

= c~ n \N+). 


Hence an h-set N is a product of two closed balls in some coordinate system 
Cat. We call numbers u(N) and s(N) unstable and stable dimensions, respec¬ 
tively. The subscript c refers to the new coordinates given by Cat- The set \N\ 
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Figure 2: An h-set N with u = 2, s = 1 and a horizontal disc in N. 

is called the support of an h-set. We often drop bars in the symbol \N\ and use 
N to denote both the h-set and its support. 

Occasionally we will say that N = {(xo,yo)}+B u (0,ri)xB s (0,r2) C R“xR s 
is an h-set. By this we will understand a ’natural’ h-set structure on N given by: 
u(N) = u, s(N) = s, c/v( x, y) = ^ x ~ xa , v ~^° ) • In the context of M 2 and u = 1, 
s = 1 we will sometimes write N = zq + [— a, o] x [—b,b\. This is compatible 
with the above convention as a defines the radius of the ball B u ( 0, a) = [—a, a] 
and b of B s (0,b) = [~b,b\. 

Definition 2 V7IT\. Def. 5] __ 

Let N be an h- set. Let b : B u (jv) —> |7V| be a continuous mapping and let 
b c = Cjv ° b. We say that b is a horizontal disc in N if there exists a homotopy 
h : [0,1] x F> u (jv) —> N c , such that 

h 0 = b c , 

hi(x) = (x, 0), WxG u(N), 

h(t , x) e N~, Wt S [0,1] and Vx £ dB u ^ N y 

The geometry of this definition is shown in Fig. [2j During the homotopy 
h deforming a horizontal disc to the map h(l,x) = (x, 0) the image of the 
boundary /i([0,1] x B s (jv)) must remain in N~. In a similar way we can define 
vertical discs. 

Definition 3 'ZCCA Def. 6] Let N be an h- set. Let b : i? s (jv) —> |AT| be a 
continuous mapping and let b c = cn ° b. We say that b is a vertical disc in N 
if there exists a homotopy h : [0,1] x B s (jv) —> N c , such that 

K = b c , 

hi(x) = (0,x), Vx 6 B s(JV ), 

h(t,x) € N^~, 'it £ [0,1] andVx £ dB s ^ N y 
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Definition 4 IZCCX Def. 7] Let TV be an h- set in M. n and b be a horizontal 
(vertical) disc in TV. We will say that x £ R n belongs to b, when b(z) = x for 
some z £ dom(&). 

By |6| we will denote the image ofb. Hence z £ |fo| iff z belongs to b. 

2.3 Cone conditions and the stable manifold theorem 

Below we present definitions and theorems that allow us to handle and ver¬ 
ify hyperbolic structures for nonautonomous ODEs using h- sets and quadratic 
forms. 

Definition 5 '/.('(", Def. 8] Let TV C R n be an h- set and Q : R ra —> R be a 
quadratic form, such that 

Q(x,y) = a(x) - /3(y), (x,y) £ R“ (A ° x R s(JV) , 

where a : —► R and /3 : R 8 ^ —> R are positive definite quadratic forms. 

The pair (TV, Q) is called an h- set with cones. 

We will often omit Q in the symbol (TV, Q) and will say that TV is an h-set 
with cones. 

Definition 6 (ZCU Def. 9] Let (TV, Q) be an h- set with cones and b : B u — > 

| TV | be a horizontal disc. 

We will say that b satisfies the cone condition (with respect to Q) iff for any 
Xi,X 2 £ B u , Xi ^ X 2 the following inequality holds: 

Q(b c (xi) - b c (x 2 )) > 0. 

The geometry of this definition is shown in Fig. [3j The entire graph of 
horizontal disc b must lie in the positive cone attached to any point z that 
belongs to the horizontal disc b (except point z). It can be shown |ZCC| that 
a horizontal disc satisfying the cone condition must be of the form 6 c (x,y) = 
(, x 7 g(x )) where g is a Lipschitz function with an explicit bound on Lipschitz 
constant coming from the quadratic form. We will use this property to control 
derivative of the parameterization of the local invariant manifold which will 
be constructed as a horizontal (vertical) disc in some h-set satisfying the cone 
condition. 

Definition 7 f ZCJCJL Def. 10] __ 

Let (TV, Q) be an h- set with cones and b : B s —» |TV| be a vertical disc. 

We will say that b satisfies the cone condition (with respect to Q) iff for any 
2/i; 2/2 G B s , yi y 2 the following inequality holds: 


Q{b c (yi) - b c {y 2 )) < 0. 













Figure 3: An h-set N and a horizontal disc in N satisfying the cone condition. 

Let us consider an nonautonomous ODE 

z' = f(t,z), (5) 

where z £ R n , / £ C' 1 (K x R”,R"). Let us denote by p(t,to,p) the solution of 
© with the initial condition z(t 0 ) = p. 

Below we present the definition of the local stable and unstable set for a 
fixed point of an non-autonomous of ODE. These are just the adaptation of the 
analogous notions for autonomous ODEs. 

Definition 8 Let zq £ R" be such that f(t, zo) = 0 for t £ R. We will call such 
zq a fixed point for @ (or f). 

Let Z C R n , zq £ Z be a fixed point for ([5| and let to G R. We define 
W t s OtZ (z 0 , p) = Wt o>z (z 0 , /) = {*: V t > 0 p(t + t 0 , t 0 , z) £ Z, mn p(t, t 0 , z) = z 0 } , 
w t 0 z( z o,<p) = Wfi z (zo,f) = \z: V t<0 </j(f + t 0) t 0 , z) £ Z, lirn <p(t,t 0 ,z ) = z 0 

u ’ u ’ — t^—oo 


When (f, f or Zq is known from the context, then we will often drop it and write 
W^ z (z 0 )orW^ z . 


Below we recall the notion of the isolating block from the Conley index 
theory with some obvious modification for nonautonomous ODEs. 


Definition 9 Let p be an autonomous local flow induced by f : R™ —> R ra . For 
S > 0 the set £ C R" is called a S— section for the flow p iff p ((—<5, <5), £) is an 
open set and the map a : £ x (—5, S) — > p((—6, S ), £) defined by <j{x, t) = p(t, x ) 
is a homeomorphism. 

For a nonautonomous ODE x' = f(t,x) inducing p we say that £ as above 
is 5-section for p for t £ T C R, iff for each to £ T the set Y, is a 5 section for 
the vector field x f(to,x) (i.e. we freeze the time). 

Let B C R" be a compact set and T C R. B is called an isolating block 
for non-autonomous ODE for t £ T iff dB = B~ U B + , where B~ and B + are 
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Figure 4: An isolating block Z for some planar ODE. The stable (vertical green 
dashed line) and unstable (horizontal red solid line) manifolds for ip inside Z are 
plotted, arrows indicate the vector field /. Dashed green and solid red border 
lines indicate the 5-sections E + and E - respectively. 


closed sets, and there exists S > 0 and two 8 — sections, E + and E such that 

B + C E + , B~ C E“, 

Mx G B + Vto G T Vt £ (—8,0) tp(t+ t 0 ,t 0 ,x) $ B, 

Vx € B~ Vto € T Vt € (0, 5 ) </?(t, x) ^ B. 

B~ and B + will called an exit set and an entrance set, respectively. 

In the present paper we will use h-sets which are isolating blocks for t > to 
(or for t < t 0 ). Simply, it means that N + and N~ are sections of the vector 
field. 


Definition 10 Let N he an h- set in IS 
for nonautonomous ODE (5) fort > t 0 , 
for t > to as in Definition 9 


n . We say that N is an isolating block 
iff and N + are 8-sections for f(t, •) 


Therefore, if N, h-set, is an isolating block, then N and N + are the exit and 
the entrance sets, respectively. 


Definition 11 Let N be an h-set, such that Cn is a diffeomorphism. For a 
vector field f on |iV| we define a vector field on N c by 

fc(z) = Dc N (cfj 1 (z))f(cf f 1 (z)). (6) 

Observe that f c is the vector field / expressed in the new variables. If / depends 
on time, then f c is time dependent too. 

Definition 12 fZCU Def. 13] 
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Let U C R" be such that U = U and int U ^ 0. Let g : U R m be a C 1 
function. We define the interval enclosure of Dg(U) by 

[Dg{U)} := |A G R nxm : G 

We say that an interval matrix [A] C R nx " (i.e. a set of matrices) is positive 
definite if for all symmetric matrices A £ [A] A is positive definite. 

The lemma below serves to motivate our definition of the cone condition for 
ODEs. 

Lemma 2 Assume that Q is a quadratic form on R n . 

Consider equation ©• Let I C R be an interval. Assume that U C R" is 
convex and the following condition is satisfied for t G / 

matrix [D z f(t,U)] T Q + Q[D z f(t,U)\ is positive definite. (7) 

Let Zi ( t ) and z 2 (t) be two different solutions of 0, such that fort G I Zi(t) G U. 
Then there exists g > 0 such that 

1 . 

j t Q(z x {t) - z 2 (t)) > 0, t£l, (8) 

j t Q{z\{t) - z 2 {t)) > ±r]Q{z]_{t) - z 2 (t)), t£l (9) 

2. if t 0 ,t + to £ I, t > 0 and Q(zi(to) — z 2 (to)) > 0, then 

Q{zi(t 0 + t)- z 2 (t 0 +t)) > e vt Q(zi(to) - z 2 (t 0 )), (10) 

3. if t 0 , t + to £ I, t < 0 and Q(z\(to) — z 2 (t 0 )) < 0, then 

Q(zi(to+t) - z 2 (t 0 +t)) < e riW Q{z 1 (t. 0 ) - z 2 (t 0 )). (11) 


. dgfix) dgfix) 
mf -,sup —- 

xeU OXj X £U OXj 


Proof: We have 

= (f(t, zi(t)) - f(t, z 2 {t))) T ■ Q ■ {zx{t) - z 2 (t)) 

+ (zi(t) - z 2 {t)) T ■ Q • (, f(t, zi{t)) - f(t,z 2 (t ))) 

= (zi(t) - z 2 (t)) T -C T ■ Q - (zi(t) - z 2 (t)) 

+ (zi(t) - z 2 (t)) T ■ Q ■ C(zi(t) - z 2 (t)) 

= ( Zi(t ) - z 2 (t)) T • ( C t Q + QC) • (zi(t) - z 2 (t)), 

where 

c = C(t,z 1 ,z 2 ) = I D z f(t,zi + s(z 2 - Zi))ds. 

Jo 


- z 2 {t)) 
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Observe that C £ [D z f(t, U)] T Q + Q[D z f(t, U)], so (8j9) follow easily. 
From © and differential inequalities we obtain 


Q(zi(t 0 +t) - z 2 (t 0 + t)) > e ±r)t Q{z 1 {t 0 ) - z 2 {t 0 )), t> 0, t 0 ,t 0 + t € I 
Q(zi(t 0 +1) - z 2 (t 0 + t)) < e ±r,t Q(zi(t 0 ) — Z 2 (to)), t < 0, t 0l t 0 + t £ I 


From the above conditions one can easily infer the remaining assertions. I 

Let us remark that from the above lemma it follows that cone condition 
(71 implies that cones C + (z) = {zi | Q(z — zz) > 0} are forward invariant for 
(51 and there is an exponential expansion in C + (z). Analogously cones C~(z) = 
{u | Q(z — u) > 0} are backward invariant and there is an exponential expansion 
in C~(z) in backward time. 

The following theorem can be easily obtained using the method as in the 
proof of Theorem 26 in [ZCCj . 

Theorem 3 Consider |5p. 

Assume that Zq £ R n is such that for any t holds f(t, zq) = 0. 

Assume that ( N,Q ) is an h- set with cones, such that Zq £ int N and which 
is an isolating block for (|5j) for t > to, cjv is a diffeomorphism and that the 
following cone condition is satisfied for t > to-' 

matrix [D z f c (t,N c )] T Q + Q[D z f c (t,N c )\ is positive definite. (12) 

Then for any t\ > to t/ie set W® N (zo,f) is a vertical disc in N satisfying 
the cone condition. Therefore there exists a function x s : {t > to} xB s -> B u , 
such that 

w t llN (zo,f) = { c f^{x s {t\,y),y)) | y £ B s }. (13) 

Analogous result is valid for the local unstable manifold. 

Observe that from the standard results about exponential dichotomies (see 
for example 0) it follows that function x s (t\,y) is smooth with respect to t\ 
and y if / and our coordinate change cn are smooth. 


3 The estimate for the derivative with respect 
to time of W^^(zo) for non-autonomous ODE 

In the sequel we represent points in R" by z — (x,y) with x £ R“ and y £ R s , 
n = u + s. 

3.1 The estimate for the time derivative of W^(zq) 

The following lemma gives computable bounds on the first derivatives of stable 
manifolds of a hyperbolic fixed point for a nonautonomous ODE with respect 
to the initial time. The proof is analogous to the proofs of the dependence 
on parameters of the stable manifolds of maps from [ iZCC| Thm. 21] with 
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some refined estimates in [WZ11 Th. 4.1] and for parameterized ODEs IISZ1 
Th. 5]. The theorem below and its proof is an adaptation of these results for 
nonautonomous ODEs. 


Lemma 4 The same assumptions as in Theorem [3| N has a natural h-set 
structure, N = B u (p x ,r x ) x B s (p y ,r y ) and that the quadratic form is given by 
(in natural coordinates) Q(x,y) = a||x|| 2 — &||y|| 2 , where a, b > 0. 

From the assumption on the cone condition for Q there exists E > 0 such 
that for t > to and z £ N there holds 


" T 1 1 Q + Q^(t,z)]v>E\\v\\ 2 , v£ 


dz 


(14) 


Put 


m = sup 

t>to,z(zN 


df, 


and 


5 = 


df, 


Of (M)J Q + Q^p(t,z) 


dt 


am 


(15) 


(16) 


Under above assumptions, if {xi,yf) £ W t s . n (zq) and ti > to for i = 1,2, 
then 

a||xi - £ 2 || 2 < S(h - t 2 ) 2 + b\\yi - y 2 1| 2 . (17) 

In particular, ifWf N (z 0 ) = {(x s (t,y),y) \ y £ B s (p y ,r y )}, then for ti,t 2 > 
to there holds 

YYl 

\\x s {ti,y) ~x s (t 2 ,y )|| < — |*i -t 2 \. (18) 


Proof: Consider the quadratic form on the extended phase-space ( t, x, y) 
Q(t, z) = Q(z) - St 2 , S > 0. 


The value of <5 will be fixed later during the proof (and it will turn out that (161 
is a good guess ). 

Observe that, if for i = 1,2 tj > to, Zi £ Wf N (f) and £,(fj + 1) for t > 0 is 
a solution of ([5| such that Zi(U) = Zi, then 


\Q{{ti+t,Zi(ti+t)) - [t 2 + t,z 2 {t 2 +t))\ < L, for t £ K+. (19) 


Indeed, since Zi(ti + t) £ N for t > 0, so \Q(zi(ti +t) — z 2 (t 2 + f))| is bounded, 
and S((t i + 1) — {t 2 +1)) 2 = 5(ti — t 2 ) 2 is constant. 

We want find a positive S, such that the positive cone defined in terms of 
Q to be forward invariant for t > to and we have an exponential growth in this 
cone, he. there exists c > 0, such that if Z\, z 2 £ N and t\ > to, t 2 > to are such 
that Q({ti,Zi) — ( t 2 ,z 2 )) > 0 and, Z\(t +1±) and z 2 (t +1 2 ) are solutions of ^ 
with Zi(ti) = Zi, i = 1, 2 , then for t > 0 holds 

Q{(ti +1, zi(h + t )) - (t 2 + t, z 2 (t 2 + t))) > e ct Q({t. i, zi) - (t 2 , z 2 )) (20) 
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if z\ (£1 + s ), z 2 (t 2 + s) G N for s G [0, t]. 

This together with (19) implies that if Zi £ W^_ N (f) and U > to for i = 1,2, 
then 

Q((h,zi) ~ {h,z 2 )) < 0 - ( 21 ) 


From this condition and the computed expression for 6 we will obtain our as¬ 
sertions. 

So now we turn to finding S. 

Let zi(t + ti) and z 2 (t + £ 2 ) be as in condition (|20|). We set 


L(t) Q((t + ti, Z\{t + ti)) — (t + t 2 , Zi(t + ^ 2 ))) — 

Q{zi{t + ti) - z 2 (t + t 2 )) - 6{h - t 2 ) 2 . 

We have 


L'( 0) — + ti) — z 2 (t + t 2 ))t=o — 

- f(t 2 , z 2 (t 2 ))) T ■ Q ■ ( 2:1 (ti) - z 2 {t 2 ))+ 

(zi(ti) - z 2 (t 2 )) T ■ Q • (/(ii, 2 i(ti)) - f(t 2 ,z 2 (t 2 ))). 

Observe that 


f{tl,Zl(ti)) - f(t 2 ,z 2 {t 2 )) = 

^((l-s)(t 1 ,zi(t 1 )) + s(t 2 ,Z2(t 2 )))ds- (ti -t 2 )+ 

j |y((l - s)(tl,^l(*l)) + s(t 2 , z 2 {t 2 )))ds ■ (Z! (ti) - z 2 (t 2 )) = 

Ct(t\ — t 2 ) + C z (zi(ti) — z 2 (t 2 )), 


where 


C t = C t (t 1 ,z 1 ,t 2 ,z 2 ) = J ^((1 - s)(ti, 2 i) + s(t 2 ,z 2 ))ds, 


C g = C z (t 1 ,z 1 ,t 2 ,z 2 ) = J ^((l-s)(t 1 ,z 1 ) + s(t 2 ,z 2 ))ds. 

From (14) we immediately obtain 

v T (CjQ + QC z )v>E\\v\\ 2 , veR n . 


From (15) it follows that 


m > \\C?Q + QC t II . 


Continuing the previous derivation we obtain 

L'(O) = (.zi - z 2 ) T (i CjQ + QC z ) { Zl - z 2 )+ 

(ti - t 2 ){CjQ)(z± - z 2 ) + (Z! - z 2 ) T (QC t ){ti - t 2 ). 
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Therefore we arrived at the following estimate 


L 1 ( 0 ) — —Q((t + t\, Z\{t + ti)) — (t + t2, Z2{t + ^ 2 )) )t—0 > 

E\\zi - z 2 \\ 2 - m|ti - t 2 \ ■ \\zi - z 2 1|. 

Now we make the use of the assumption — {t. 2 , Z 2 )) > 0. We have 

^ll~i - z 2 \\ > ||xi - x 2 || > |ti -1 2 |- 
Therefore we obtain 


l\ 0) > e\\z 1 - 2t 2 || 2 - ll-i - z 2\\ 2 = [e - 


z 1 - z 2 


Hence the positive cone for Q will be forward invariant provided 


which implies that 


E > m\ 

V 0 


s> 


E 2 


( 22 ) 


Let us fix any 6 satisfying (22). 

From the above reasoning we know that the positive cone for Q will be 
forward invariant and if Q{{t±,Zi) — (^2,^2)) > 0, then for any t , such that 
Zi{U +t) £ N holds 

^Q((f + t 1 ,z 1 (t + h)) - {t + t 2 ,z 2 (t + t 2 ))) > ^E - m^J^j \\zi - z 2 \\ 2 > 

cQ((ti + t, Z\(t\ + t )) — (^2 + t, 22(^2 + 1 )) 

for some positive c. Therefore 

Q((t + ti,zx(t + h)) - {t + t 2 , z 2 {t + t 2 ))) > Q((h,zi) - (t 2 ,z 2 ))e ct 
for t > 0 such that both z\{s + ii) and Z2(s + t2) are in N for s £ [0, t]. 


Therefore we have established (21) for any S satisfying inequality (22), but 
passing to <5 —>• we obtain it also for 5 = g ^-. 

If we set Zi = ( Xi , yi), then from (21) we have the following inequality 

a\\x\ - x 2 \\ 2 < 5(ti - t 2 ) 2 + b\\yi - y 2 \\ 2 . 

This proves Inequality is obtained, when we apply fl7| ) to Zi = 

( ti,x s (ti,y),y). 
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In the application of Lemma [4] in this paper we will only use estimate (18). 
Moreover, since we are using the natural coordinates on N cone condition (12) 
becomes 

matrix [D z f(t,N)] T Q + Q[D z f(t,N)\ is positive definite. (23) 
This simply means that we do not need to use the coordinate change c/v. 


3.2 Estimate for time derivative of W^ N (zo) 

By reversing the direction of time from Lemma [I] we obtain the following state¬ 
ment about the derivatives of the unstable manifold. 

Lemma 5 Assume that f(t,zo) = 0 for all t £ R. Let N has a natural h-set 
structure, N = B u {p x ,r x ) x B s (p y ,r y ), Zq £ intlV and that quadratic form Q 
on is given by (in natural coordinates) Q(x,y) = a||cc|| 2 — 6||y|| 2 , where a, b > 0. 

Assume that for t < to the set N is an isolating block for 0 and f satisfies 
cone condition on N with respect to quadratic form Q. 

From the cone condition for Q there exists E > 0 such that for t < to and 
z £ N 

yT + v - -^IMI 2 ’ ' ,eK " 

Put 

m = sup 

t<to,z£N 

Under above assumptions, if Wf N (zo) = {(x,y u {t, x)) \ x £ B u (0, 1)}, then 
for t\,t -2 < to holds 

TYl 

\\y u {h,x) - y u (t 2 ,x) || < — |ti - t 2 |. 

Proof: Transformation t —► — t changes Q —► — Q and (x,y) —► {y, x). The 
quantities in Lemma E] are not affected by this change. I 


df , 


Q + 


4 Estimates for the stable manifold for r 
for equation ([l]) 

In the limit r —> oo equation (jTJ) becomes 

u” = u — u 3 . 


OO 


We consider the first order system corresponding to (241. 


v! — V, v' = u — u 3 . 


(24) 

(25) 


It is easy to see that (25) is Hamiltonian with H{v,u) = \\ 
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It is clear that (0,0) is a hyperbolic fixed for ( [25] ). Its stable and unstable 
manifolds coincide to form an eight-shaped loop (called a separatrix) shown in 
two upper row drawings in Fig. [T] Any solution of ([Tj) which converges to 0 for 
r —► oo when plotted on (A, A') plane should approach this separatrix. In Fig. [l] 
this is shown for six orbits. 

We introduce coordinates, which diagonalize the linear part of (25) by 


x e = u + v, 
y e =u-v. 


(26) 


The inverse is given by 


u = (x e + y e )/ 2, 
v = (x e - y e )/2. 


In new coordinates system (25) becomes 


x’ e = x e + g(x e ,y e ), 
y' e = z/e - g(x e ,y e ) 


(27) 


where 


. . x — y x + y 1. ,o 

g{x,y) = -—:- + —w- - -(x + y ) 3 . 


2 r 


8 r 2 


4.1 Isolating block 


We look for the isolating block N e = [—d \, df] x [—d|, d|] f° r (27). 
Lemma 6 Put d = d\+ d\ and assume that 


min{d®, d|} > d 


1 _ d 2 

2r + 8r^ + IT 


(28) 


Then N e is an isolating block for (21). 


Proof: 

First we verify the exit condition. For |x e | = d\ we have 


x e x e X e T x e g(x e , y e ) 


2 1 x e -y e (x e + y e )x e (x e + y e ) 3 x e 

Xe 2 r Xe ^ 8r*-8- 

2 1 d\d d\d d 3 d\ 


> 


^-2-r-8^- 


d\\d^d [- + - 2 


<p_ 

8 


Therefore x e x' e > 0 when (28) is satisfied. 
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In a similar way we verify the entry condition. For \y e \ = we have 


- (v'eVe) 


y 2 e + y e g(x e ,y e ) 


> 



1 

8 r 2 



Therefore — (y e y e )' > 0 when (28) is satisfied. 


4.2 Cone condition 


Our cones are defined in terms of the quadratic form 


Q 


a 0 
0 -1 


(29) 


Let 

C(NJ) = [Df T (N)]Q + Q\Df(N)]. 

where [Df(N)] is defined as follows: M £ [Df(N)] iff there exists a pair of 
points Zi,Z 2 £ N (can be equal) such that M = f Q Df{z\ + t(z 2 — zi))dt. 

Lemma 7 Assume that N is a convex and compact set and f : is C 1 

smooth. Assume that for every z £ N there exists e(z) > ejy > 0, such that for 
all v £ R" there holds 


v T (Df T (z)Q + QDf{z))v > e(z)||w|| 2 . 
Then for every M £ C ( N , /) and v £ R" there holds 
v t (M t Q + QM)v > eAr||u|| 2 . 


Proof: For z\,Z 2 £ N and v £ R n we have 


J Df T (zi + t(z 2 - zi))tdtj Q + Q (^J Df(zi+t(z 2 -zi))dt 

[ v T (Df T {z l + t(z 2 - zi))Q + QDf(zi + t{z 2 - zf)))vdt > e N 
Jo 


V = 


Lemma 8 Let N e = [—df,df] x [— , cZf] an( ^ Q be given by (29) with a = 1. 
Assume that for some E > 0 and r* > 0 there holds 

1 > — + -(d® + dt,) 2 + E/2. 


Then for all v £ R n , M £ C(N e , f) and r > r* 

v t (MQ + QM)v > -LIMI 2 , 


(30) 


i.e. the cone condition holds on N e for r > r». 
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Proof: According to Lemma[7]it is sufficient to show that the matrix Df T (x e , y e )Q+ 
QDf(x e ,y e ) — E ■ Id is positive definite for each (x e ,y e ) £ N e . Put 


Mi = 
We have 


1 0 
0 -1 


, M 2 = 


1 1 
-1 -1 


1 3 


, m 3 = 


-1 1 
1 -1 


Df{x e , 2/e) = Mi + ( - - (x e + 2/e)“ ) M 2 + — M 3 . 


Easy calculations give 

M?Q=\“ ? 1 , MjQ = 


a 1 
a 1 


, MlQ = 


2 r 


—a —1 
a 1 


We obtain 

Df T (x e , y e )Q + QDf( x e , y e ) = 


2 a 0 
0 2 


+ ( ^ ~ | (*■ + v^ 


2a 1 + a 
a A 1 2 

1 

2r 


—2a a — 1 
a- 1 2 


Let us denote 

W = 8r2 f^e + ye) 2 - (31) 

Then 


Df T (x e , y e )Q + QDf(x e , y e ) - E ■ Id = 

2a(l + tu--E u>(a+l) + ^=^ 
w(a + 1) + ^ 2r ^ 2 (l + W + — E 

It is immediate that for the positive definiteness of this matrix we need 

1 

1 + re — --£ > 0. 

2 r 

From now on we set a = 1. The Gersgorin Theorem gives us the following 
sufficient condition for the positive definiteness of the matrix Df T (x e ,y e )Q + 
QDf(x e ,y e ) -E-Id 


2 

2 


1 + re — 


2 r 


1+ " + 2; 


— E > 2|io|, 

— £ > 2\w\. 
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We end up with the following condition 


1 + w ~ 7T > M + E/2. 
Zr 


Substituting back for w defined by (31) we obtain two inequalities that must be 
simultaneously satisfied 


1+ ^2 “ g^e+Z/e) 2 - Y r > ^2 + g^ Xe + y^ 2 + E / 2 ’ 

1-— > ~(x e + y e ) 2 + E/2. 


Since for (x e . y e ) £ N e there holds \x e + y e | < (df + d|) we obtain the following 
condition 

1 > — + - (d® + d^) 2 + E/2 


and Lemma [T] implies the assertion (301. 


4.3 Estimate on constant m 

Lemma 9 Let N e = [—df,df] x [—d|,d|] and put d = df + d|. Then for 
z = (x e ,y e ) £ N e and r > r* there holds 


Q + Q?f(r,z) 


^( 1+ J 


Proof: Since we assumed that Q is diagonal it is easy to see that 


Q + Q%( r ’ z ) 


= 2 


z ) 


We have 


hence 


df, , f x e ~y e 3 




i 

- 1 


Finally we obtain 





























4.4 Some numbers 


Consider equation ( |27[ ) and quadratic form ( |29| with a = 1 (this is the value 
used later in the actual proof). Lemma [6] and Lemm a [8] guarantee that the 
set N e = [—d®,df] x [—d|,djj] is an isolating block for (27) satisfying the cone 
condition for all r > r* provided the following inequalities hold true: 


min{d®,d|} > 


,1 1 d 2 

d| ^: + 8 ^ + ¥ 


£ 1 3 , 2 

2 + 2^ + l d 


< 1 , 


where d = df + d§. 

For example, if r* = 2 and d\ = d| = 1/4 then the first inequality holds 
true. From Lemma 9 we get the bound m < Now the constant E should 

satisfy E < 2 — 3/2cP — 1/r* = 9/8. Take for instance E = 1. This gives us a 
bound for the derivative |x'(r)| < ^ ~ 0.121534. 

For r large the second condition gives an upper bound for d < | \/3. 


5 Estimates for the unstable manifold for p —> 
—oo for equation ([l]) 

Let us rewrite here system Q 

w" — w/A = e 2p (w — w 3 ). (32) 


We want to study the unstable manifold of (0,0) for p < 0. 

Passing formally to the limit p —> —oo in (32) we obtain the following linear 
system 


w” - w/A = 0. 

Let p = w' and let us introduce new coordinates Xb, yh 

w = x b + y b , 
x b - Vb 


The inverse is given by 


P = 


x h = 


Vb = 


2 ' 

w + 2p 
2 ’ 
w — 2p 


In the coordinates (xb,y b ) (32) becomes 


1 


x b = 2 Xh + e2p d(x b,Vb), 
Vb = ~2 Vh ~ e2p 9( x b,Vb), 


(33) 


(34) 


(35) 
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where 


g{x h ,y h ) = (x h +y h ) - ( x h + y h f. 


5.1 Isolating block 


Lemma 10 Let N b 

(35) for any p. 


[—cZ, cZ] 2 . If d < 1/2 then N h is an isolating block for 


Proof: Observe that if |xb + ?/b| < 1 then g(xb,yb) has the same sign as x^ + yb- 
The sign of Xb+2/b coincides with the sign of the dominant coordinate. Therefore 


x h g(x h ,y h ) > 0 for |cc b | = d, 
ybg(x h ,y h ) > 0 for \y h \ = d. 


This implies that 


x'bXh = T^b/2 + e 2p x h g(x h ,y h ) > ^b/ 2 if |*b| = d, 
-VbVb = \yll% + e 2p y h g{x hl y h ) > \yl/2 if |y b | = d. 


5.2 Cone condition 


In this section we will derive inequalities that guarantee the cone condition on 


an isolating block for (35). Let us denote by / the vector field (35). We have 


Df{xb,yb ) = 


0 


1 
2 

o , 


+ e 2p (l — 3(a:b + y b ) 2 ) 


1 1 
-1 -1 


Put 


Q = 


a 0 

0 -1 


(36) 


Lemma 11 Assume that the isolating block for (35) has the form N h = [— di,d\]x 
\-d%,d\]. Put 


w = e 2p 


(1 - 3(x h + y h ) )■ 


If 


1 + 2 w > E> 0, (37) 

a(l + 2w) 2 - E(a+l){l + 2w)+E 2 > w 2 (a + 1) 2 (38) 

for {xb,yb) G N h and p < po then the cone condition is satisfied on N h for all 
p < po for \35ty with the constant E. 
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Proof: Denote by / the vector field (351. Direct computation gives 


Df T {x b , y h )Q + QDf{x b , 2/b) = 

a 0 
0 1 

a(l + 2 w) w(a + 1) 
w(a + 1) 1 + 2 w 


+ e ip (l-3(x b +2/b) ) 


2a a + 1 
a + 1 2 


We will apply the standard criterion for positive definiteness of matrices to the 
matrix M = (Df T (xb,yb)Q + QDf(xb,yb ) — E ■ Id). We want M 22 > 0 and 
det M > 0. The first condition gives ( [37| , the second reads 

(a(l + 2 w) — E) (1 + 2 w — E) — w 2 (a + l) 2 > 0 


which is equivalent to (381. 1 

Using the Gersgorin theorem mm we can derive conditions, which might 
more suitable for easy estimates to be done by hand. The positive definiteness 
of 

a(l + 2u>) w{a + 1) 
u>(a+l) 1 + 2 w 

is implied by the following two inequalities 


a(l + 2 w) — w(a + 1) > E, 
1 + 2w — w(a + 1) > E. 


(39) 


This leads to the following lemma. 

Lemma 12 Assume that (d^ + d^) 2 < 0 < a < 1 and the following condition 

is satisfied 


i(e 2po + 1) - e 2po > E. 


Then the cone condition is satisfied on N h for (35) for p < po with the constant 

E. 


Proof: Under the assumption on d\ and df( we have 

w = e 2p (l - 3(a: b + 2/b) 2 ) > 0. 

Since 0 < a < 1, the positive definiteness of Df(xb,yb) T Q + QDf(xb,yb) is 
implied by only one condition (39) which reduces to 

a(w + 1) — w > E. 


Observe that the range of w for (xb, 2/b) £ N h and p < p 0 is given by 0 < w < 
e 2po = wq. Since a < 1, the function w 1 —> a(l + w) — w is non-increasing, hence 
it is enough to have this inequality for the largest possible value wq = e 2po . 


23 














5.3 Estimation of constant m when p —> —oo. 


Le mma 13 Denote by f the vector field (35). Assume that the isolating block 
for (35) has the form N h = [—d\,di] x [— dfficfi)] with d := max{< 1/2. 
Assume the quadratic form is diagonal as in (36). Then for z = (xb,yb) £ N h 
and p < po there holds 


Q+Q %^ z) 


< 8 de 2po \Ja 2 + 1. 


Proof: For 2 = ( Xb,yb ) £ N h we have 

9f, 


dp 

Since Q is diagonal we get 


0, z) = 2e 2p g(x h ,y h ) 


1 

-1 


df 

= 2e 2p g{x h ,y h ) 


Given that 

\g(x h ,y h )\ = |a.’ b +j/b|(l- (a"b + 2/b) 2 ) < \xb + Vb\ < 2 d. 
for (xb,yb) £ N h we obtain 

Q (p, z) < Ade 2po \]a 2 + 1. 
dp 


5.4 Some numbers 


We will show now that we can construct a quite large isolating block on which 
the cone condition is satisfied and still we can have good estimates on the 
derivative on par am etrization of unstable set with respect to p. 

the set N h = [—d, d\ 2 is an isolating block provided d < 1/2. 
guarantees the cone condition on N h with constant E for p < po 


From Lemma 
Lemma 


12 


provided the following two inequalities are satisfied 


d 


1 


(e 2po +1) - 


O 2po 


2 ^ 3 ’ 
> E. 


In the case a = 1 (this is the value used in our computer assisted proof) the last 
inequali ty r educes to E < 1 and po can be an arbitrary number. Finally, from 


Lemma 13 we obtain an estimate on the constant m = 8de 2po V2. 


Explicit numbers satisfying all required inequalities: 
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• Put d = 1/4, po = In 2, a = 1 and E = 1/2. Then we have m = 8a/2, 
Wu(p)\ < m/E = 16\/2. 

• Put d = 1/4, po = ln2 -4 , a = l and E = l/\/2. Then we have m = 
V2/128, \y' u {p)\ < m/E = 1/64. 


6 Proof of Theorem |T| 


In this section we give a computer assisted proof of the existence of six geomet¬ 
rically different connecting orbits for equation ([l]), i.e. solutions satisfying ([2]). 
Each orbit makes different number of (half) revolutions around the separatrix of 
equation ( [25] ) before approaching the equilibrium point {A, A') = (0,0). Three 
of them converge to (0, 0) from the right side {A > 0) and three of them from 
the left side {A < 0) - see Fig. [l] and Fi g. [5| 

Before we give the proof of Theorem [ljlet us make several remarks about its 
content. 


• W u can be expressed using other independent variable p = In r and w(p) = 
A{e p ) (see Section [HJ) as follows. 


W u = {(po, wq, w' 0 )) £ R x ffi. x 


such that the backward solution of (32) w(p) with initial condition 
w(p 0 ) = w 0 , w'(po) = w o satisfies lim p _ > _ 00 (u;(p), w'{p)) = 0}. 


In fact in the proof we use the above description of W u . The expression 
for W u given in the assertion of the theorem is obtained after we return 
to the original variables r and function A(r). 

• The statement about the number of local extremes means that the trajec¬ 
tory r —»• (A n (r), A' n (r)) intersects transversally the axis A' = 0 exactly n 
times. 

In Fig. [5] the six connecting orbits are shown in two time scales: r for r > 1 
and p = lnr for r € (0,1). For r = 1 they coincide, i.e. 

{M r = 1 )M\r = 1)) = {A{p = 0), A'{p = 0)) 


therefore they appear as continuous curves. This figure illustrates the conver¬ 
gence of orbits to (0,0) when p —► — oo. 

The proof of Theorem [l] is based on the theorems introduced in the previous 
sections. Using a non-rigorous simulation (the bisection) we found good initial 
conditions for the connecting orbits. Then using validated numerics we verified 
the existence of isolating blocks, cone conditions and conditions for the inter¬ 
section of the stable and unstable manifolds. The details will be given in the 
next sections. 
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Figure 5: The six transverse connecting orbits resulting from Theorem [I] The 
trajectories are plotted in two time scales. For r > 1 the points (A n (r), A' n (r)) 
are shown while for r £ (0,1) we plot ( A n (p ), A' n (p)) with p = lnr. This figure 
illustrates the convergence of orbits to (0,0) when p —► — oo. The figure-eight 
curve in the middle is the separatrix for (25). 
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In Sections [4] and [5] we have found analytically the isolating blocks on which 
the cone conditions are satisfied for the begin and end our connection. How¬ 
ever, using validated numerics we can find much tighter bounds for the unsta¬ 
ble/stable manifolds and similarly better bounds for the time-derivative of the 
parametrization of unstable/stable sets. 

The input arguments to the algorithms that verify the existence of connecting 
orbit are the following 

• n - number of intersections of the trajectory with the line A! = 0, 

• r n - an initial time, 


• A r n - time range around f n . 


These parameters are listed in (41) and found in nonrigorous simulation (bisec¬ 
tion). 


6.1 Isolating block and cone condition for p -A — 00 . 

By (*bi2/b) we denote coordinates that linearize the system at p —> —00 as 
defined in p4| . We define an isolating block 

N h = [~d\,d\] x \-d\,d\] = [-0.125,0.125] x [-2.8 • KT 8 ,2.8 • 10" 8 ]. (40) 


Let us define the following constants 


n 

fn 

A r n 

1 

0.003288250 

4-10" 7 

2 

0.001184020 

6•10 -8 

3 

0.000650050 

3•10 -8 

4 

0.000424204 

2•10" 8 

5 

0.000304427 

1•10" 8 

6 

0.000232050 

8•10" 9 


Lemma 14 Put p * = In (f 1 + Ari), where fy and Ar! are defined in (fl). 

• For p < p* set N h is an isolating block for the system (35). 

• For p < p* the cone condition is satisfied on N h with quadratic form 

Q{x b, 2/b) = — y\ and thus for all p < p* the set 


W p,N h (°i °) = {(xb,y u (p,Xb)) \x h e [-^1,^1]} 


is a horizontal disc satisfying the cone condition. 


• The following estimate holds: 


dy u 

dp 


(p,x b) 


< 7.65 • 10" 6 . 
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Proof: Let f(p,x b ,yb) = (fi(p, x b , y h ), f 2 (p, x h , y h )) be the vector field as 
defined in (35) and put f p = f(p, Let d\,d 2 be the sizes of N h - see (40). 


Direct evaluation in interval arithmetics gives the following inequalities 


/l((-0O,p*] 

X {4} X [- 

-d h 2 ,d h 2 ]) 

> 

0.0625 > 0, 

/i(( oo,p*] X 

{-<*?} X [- 

AA\) 

< 

-0.0625 < 0, 

/ 2 ((-oo,p»] 

X [-d?,4] 

X A}) 

< 

-4.8127929800192025 • 10~ 8 < 0 

/ 2 ((-oo,p»] X 

[-d?,d?] X 

{-d h 2 }) 

> 

4.8127929800192025 • 10” 8 > 0. 


This proves that N h is an isolating block for all p < p*. 

Let us fix p £ (—oo, p*] and recall that the quadratic form we are using is 

Q(x b, yh) = x b^ Ub- Since d\ + d 2 < A from Lemma 
condition is satisfied with any constant E < 1. 
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we get that the cone 


Finally, from Lemma 13 and Lemma[5]and passing to the limit with E —> 1 
we get that 


dy u , , 

— P ( ftIb) 


< 


< 


,P*\ x N h ) T Q + <3§{((-oo,p*] x N h ) 


dp ' 


E 


AV2e^\d h 1+ d h 2 \ E-+1- ]d b + d bj 


= 540761062255450812 • lO" 23 ^ < 7.65 • 10" 6 . 


6.2 Isolating block and cone condition for r —>■ oo. 

Let us consider an isolating block for r —> oo. Put 

N e = [-di,di] x [-(%,<%] = [-0.0015,0.0015] x [-0.01,0.01], (42) 


Lemma 15 Let r* = 6. 


For r > r* the set N e is an isolating block for the system (25). 


• For r > r* the cone condition is satisfied on N e with quadratic form 
Q(x e , y e ) = x\ — y\ and thus for all r > r* the set 


W'jv.(0,0) = {{x s (r,y e ),y e ) \y e £ [-d 2 ,d 2 ]} 


is a vertical disc satisfying the cone condition. 


• The following estimate holds: 


dx s 

dr 


(P, 2/e) 


< 0.000252. 
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Proof: Let f(r,x e ,y e ) = (fi(r,x e ,y e ), f 2 (r,x e ,y e )) be the vector field as 
defined in (25) and let d \, d% be the sizes of N e - see (42). Direct evaluation in 
interval arithmetics gives the following inequalities 


/i([r*,oo) x {c^} x [~d e 2 ,d e 2 \) > 

/i([r*,oo) x {-d®} x [-d e 2 ,d e 2 \) < 

/ 2 ([r*,oo) x x {d®}) < 

/ 2 ([r*,oo) x [~d®,d®] x {-d 2 }) > 


0.00051215277777777761 > 0, 
-0.00051215277777777761 < 0, 
-0.010737706706597221 < 0, 
0.010737706706597221 > 0. 


This proves that N e is an isolating block for all r > r*. 

From the proof of Lemma [8] it follows that the cone condition is satisfied on 
the set N e for r > r* and with a = 1, E > 0 provided the following inequality 
holds true 

E <2- — - hdl+dl) 2 


Substituting r* = 6 and d \, d 2 form (42) we get the upper bound 


E < En — 


43995239 


1.8331349583333332. 


24000000 

Finally, from Lemma [9] and Lemma [5] we get that 


dx s 

dr 


{r,y e ) 


([r*,oo) x N e ) T Q + Q§y([r*,oo) x N e ) 
E 


< 


V2(dt + d§) 

r 2 
' * 



Passing to the limit with E ^ E 0 we obtain 


dx s 

dr 


(r, Ve) 


< 


V2{d\ + d| 


1 + 5 ? /E " 


140875 

791914302 


V2 < 0.000252. 


6.3 Shooting between manifolds. 


Let us define three Poincare sections. The first section is expressed in coordi¬ 
nates (r,x b,2/b) - see (34) - that linearize the system at p —> — oo. The section 
contains one of the exit edges of the set N h , namely 


II b = {(r,x h ,y h ) | x h = d b } . 

We will use (r, y^) coordinates to describe points on II b . 

The two remaining sections are expressed in coordinates that define the set 
N e and each of them contains one of the entrance edges of this set 


II± = {(r,x e ,y e ) | y e = ±d%} ■ 
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We will use (r, x e ) coordinates to describe points on IIj_ when the sign will be 
clear from the context. By P_ : II b —> III and P + : II b — > 11+ we denote two 
Poincare maps. By P±(r, y b ) we denote the n-th itersection of the trajectory of 
|lj) starting at (r, d\,y b ) with the section ITj_. In what follows we will always 
use odd number of intersections for the mapping PI and even for the mapping 
P". Therefore we will skip the subscript ± to simplify the notation and write 
just P". 


Lemma 16 Let r n , A r n , n = 1,..., 6 be as defined in (41 ) and put 


r~ = r n - A r n 
r]\ = r n + A r n 

P 2 = [ — ^21^2]' 


Then for n = 1,..., 6 the mapping P” +1 is well defined and smooth on the set 


K, r n] x [-d 2 ,d 2 ]- 


Moreover 


pn+1 ({ r n} x D 2 ) c {(r,x e ) : |x e | > d\,r > r* = 6}, 
n x P n+1 (r~,y h ) ■ -jr x P n+1 (r+,y h ) < 0 for y h G P b 

and each trajectory of a point (r^, j/b) € [r“,r+] x [— intersects the axis 

A! = 0 exactly n times before reaching the section II?-^. 

Proof: The proof is computer assisted and uses an rigorous ODE solver from 


the CAPD library jCAPD| . In (43) we give rigorous bounds for the return time, 
i.e. for all (r b ,2/b) £ [?A,r+] x [-d^,d\] the image (r e ,x e ) := P (n+1) (r b , y b ) 
exists and r e belongs to the corresponding interval listed in (43). This shows 
that such an intersection occurs always for r e > r* = 6. 


n 


bound for the return time r e 


[6.5694270711914049. 
[9.547364685097655, 
[12.63188037430908, 
[15.467339314615232. 
[18.659449295387446. 
[21.645791630860828. 


, 6.8663028711914071] 
9.8754898050976578] 
12.975630434309084] 

( 15.811089354615238] 
, 19.009330830371958] 
,21.989541646860836] 


(43) 


In (44) we list estimates which establish the remaining inequalities. We see that 
the ^-coordinate of P n+1 ({r±} x D b ) h as opposite sign for r~ and rif and in 
each case its absolute value is bigger than d\ = 0.0015. 


n 

7r*P" +i ({r-} x D b ) 

7r x P"+ i ({r+} x D b ) 

1 

—0.0032[19,26] 

0.0032[89,97] 

2 

0.00189[2,8] 

—0.0018[86,93] 

3 

-0.0019[87,92] 

0.0019[37,43] 

4 

0.00209[5,9] 

—0.00221 [0,5] 

5 

—0.00155[5,9] 

0.0015[78,82] 

6 

0.00173[1,4] 

—0.0016[57,61] 
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- 0.03 - 0.02 - 0.01 0.00 0.01 0.02 0.03 


Figure 6: The set N e , Poincare section 11+ and an indication that P^(r 1 ,D\) 
and Pl(rt , D\) are mapped to opposite sides of the set N e . 


In Fig. [6] we show the isolating block N e , Poincare section 11+ and by black dots 
we marked the image P+({rf} x D b ). i 


Proof of Theorem [TJ 

Recall that we used (r b ,y b ) coordinates on section II b and ( r e ,x e ) coordi¬ 
nates on 11+ . From Lemmas [141 and [l5l we know that 


W P >(0,0) = {(x b ,y u {p,x b)) |x b G 
IF r s iJV e(0,0) = {(x s (r,y e ),y e ) \ y e G [—<i|, d§]}, 


where y u (-,'), x s (-, •) are smooth functions. Put 


y u (rb,x b ) := y u (lnr b ,x h ) 


and consider the mapping 


F n (r h ) =t T Xe P n+1 (r b , y u (r b , d b )) - x s (-K r P n+l (r h , y u (r b , d b )), d^). 

Zeros of this function correspond to connecting orbits we are searching for. 
We will show that for n = 1,..., 6 this function has a unique zero in [r~ 

is continuous and F n (r. 


From Lemma 16 


:) F n{r. 


n 5 
+ 


r 1 
' n , 

) < 


o, 


the function F n 

because \x a (n r P n ~(r b , y u (r b , d b )), d|) | < df. This proves that F n has zero in 
[ r nPn]- For uniqueness it is enough to show that F^(r) ^ 0 for r G [r“,r+]. 

The proof is computer assisted and we will give details for n = 1, only. For 
n = 2,... ,6 we will give the necessary estimates. Using the CAPD library 
( API)" we computed a rigorous bound for the derivative of the return map P+ 
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on the set W = [r 1 ,rf ] x [—d^d^] and we got 


DPl(W) = 


(WO 


dr T P•* 
dr b 

dn *°E±(W) 


dir r P_ 


D 2 

HW) 


c 


dr b \ rr > dy b 

-292366,565586] [-58.3544,88.8462] 

[1538.59,16826.7] [-0.253485,2.36953] 


From Lemma 14 and Lemma [15] we have bounds for the partial derivatives 
of the parametrization of invariant sets with respect to time variable, namely 

< 7.65 • 10 -6 , 

0.000252. 


9y u , , b . 


Taking into account time rescaling p = In rb we obtain the estimate 


d 

dr h 


y u {r h ,d\) 


d 

dr h 


y„(lnr b ,d0 


< 


7.65 • 10 


-6 


51 

21919 


< 0.00233. 


Put [x' s ] := 0.000252 • [-1,1] and [y' u ] := 7.65 • 10" 6 • [-1,1], 

Gathering this together we obtain a rigorous bound for the derivative 


]) C 


^ £± (W) + ^^(W)[y' u : 


- x‘ 


C 


d ^ p +(W) 


dn r Pl 


mm) 


c 


dr b 

•s] ( drh V > 9y b 

[1538.59,16826.7] + [-0.253485,2.36953] [y' u 
[a/J ([-292366,565586] + [-58.3544,88.8462] [$J) 
[1396,16970], 


(45) 


which is nonzero (see Remark 17 for short comments about the possibility for 
obtaining more precise bounds). 

In a similar way we computed this derivative for n = 2,..., 6 and we got the 
following bounds 

[-56203.7,-12850.2], 

[16342.1,137481], 

[-183339,-24681.3], (46) 

[39626.1,307922], 

[-435119,-1138.93]. 

It remains to show the transversality of A n . 

Observe that till now we established the following facts 

• Let W^ p<p j Nb C (—oo, p*] x N h be a set consisting of all points (po, • 2 o) ; 
such that there exists a solution z(p) for p < po of equation (32) (in 


mr-2 


c 

mrs 

,^3 + ]) 

c 

mu 

,tf]) 

c 

mu 

) r 5"]) 

c 

mu 

,rf}) 

c 
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coordinates in which N h is defined this is the system (35)), such that 
z(po) = zo, z((-oo,p 0 ] C N b and linip^-oo z{p) = (0,0). 

We have proved that there exists a smooth function y u (p,x b), such that 


W (p<p*),n b = {( P,Xb,Vu{p,x h )), | p G (— 00 , p*\,Xb G [-d b ,d b ]} 


Let W( r>r ^ NB C [r*, 00 ) x N e be a set consisting of all points (ro, £ 0 ), such 
that there exists a solution z(r) for r > ro of equation (25) (in coordinates 
in which N e is defined this is the system ©), such that z(r 0 ) = zq , 
2 [ro,oc) C N e and lim r _*x> z(r) = (0,0). 

We have proved that there exists a smooth function x s (r, y e ), such that 
W (r>r,),N* = {(r,x s (r,y e ),y e ), \ r G [r*,oo),y e G [~d e 2 ,d e 2 ]}. 


The sets W™ K ^ Nh and W^ r>r ^ NB are two dimensional surfaces in the ex¬ 
tended phase space RxR 2 (or [ 0 , 00 ) x R 2 depending whether we use the vari¬ 
able p or r). These surfaces can be ’globalized’ to invariant sets W u and W s , 
by applying to W^ p<pf) Nb and W^ r ^ Ne the flow induced by Q. In this way 
we obtain two-dimensional immersed invariant manifolds. 

Each of the solutions A n whose existence we already established belongs 
to W u D W s . We want to prove that the intersection along such solution is 
transversal. 

Let 7 ( 7 ") := (r, A n (r), A' n (r)) for some n = 1,..., 6 be one of our connecting 
orbits. Let pb = (rt>, d\ , y u (fb, d\)) and p e = (f e , x s (r e , d|), d|) be the points on 
7 on sections II b (i.e. Xb = d b ) and II e (i.e. y e = d|), respectively. To prove 
the transversality it is enough to show that the tangent spaces to W u and W s 
at p e satisfy 

T P W U + T P W S =R 3 . 

Since the point p e is in the isolating block N e in coordinates ( x e ,y e ) we have 

T p W s = T Pe W[ r > rt)iN . = span { 7 '(f e ), (l, ^(f e , d£), 0 ) | 

and these vectors are linearly independent because the last coordinate of 7 '(f e ) 
is nonzero. In the above formula the vector 7 '(f e ) is the direction of the vector 
field (and the curve 7 ). The second vector 

91 = ^ 1 , 

spans T Po (w^ r>r ^ Ne HIT^, which is one dimensional. 

T p W u is obtained from T Pb W^ p , ) Nh , as follows. T Pb W { p < p -),N h is s P anned 
by 7 , (n>) and 

92 = (l, 0 , ^ r {fb,d\)^j , 
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which spans the one-dimensional space T Ph [W™ p<p ^ Nh D II b j. 

Put t = f e — fb which is the transition time for P ra+1 at p e . The shift along 
the trajectory of 7 by r maps "/'(fb) to 7 '(f e ). The Poincare map P” +1 acts on 
g 2 to produce a vector <72 as follows 


= 


7r a e 52 = 


%a52 = 

We have 


pn+1 a pn +1 

g r (fb,y u (f b ,d\)) + ^ (f b; ^(r- b ,d b ))^(f b; ^), 

~ ,b\\ , STr^P”-^ 1 _ b \\dyut- ib\ 

- Q~ r - (r h ,y u (r b ,d()) + —— - (r h ,y u (r h ,d^)) — (r h ,d^), 

0 . 


T Pb W u = span {"/’(r e ), 32 } ■ 


Therefore we need to show that the vectors 7'(f e ), gi, <72 are linearly indepen¬ 
dent. Observe that 7r y B "/'(f e ) ^ 0, because the section II e is transversal, and 
7r yo5i = ^vJ)^ = 0- Therefore, it is enough to check that the projections of 
< 7 i ,<72 on (r,x e ) subspace are linearly independent. This means the following 
determinant is nonzero (we dropped arguments in the partial derivatives) 


D = det 


1 dn ,.P ,l+1 1 dir r P n + 1 dy n 

5 dr ' dy dr 

dx B d^ e P n+1 d^P" +1 dy u 

dr ’ dr dy dr 


Observe that this leads to the following condition 


D = 


dn XB P n+1 , dn XB P n+1 dy u 


dr 


dy 


dr 


dx e (dn r P n+1 
dr \ dr 


dn r P n+1 dy u \ 
dy dr J * 


We have 


dn x P n+1 
dr 

dn r P n+1 

dr 


dTT Xe P n+1 dy v 
dy dr 
dTT r P n+1 dy u 
dy dr 


= —p7r x P n+1 

dr 

- — — 7 r P n+1 

— n r r 

dr 


(r,y u (r,d\)), 
(r, y u (r,d\)). 


therefore 

D = F' n (f e ) ^ 0 

hence the intersection is transversal. 1 


Remark 17 The estimation of F[ given by (45) is very rough, but it is sufficient 
for our purposes. It is easy to see that the main reason for such big diameter 

This number can be 


dr h 


of our estimate for F[ comes from our bound for 
non-rigorously estimated by taking finite differences using data from (44) (first 
row) and Ary from (41) to obtain 


TT Xc P 2 (r+,D%) - tt Xb P 2 (x 1 ,£%) 
2Ary 


(0.003294 + 0.003223)/(8 • 1CT 7 ) = 8146.25. 
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Subdividing Ar± onto 100 pieces we got a rigorous estimate on the derivative 
F{([rr,r+]) C [7902.54,8399.24], 

From the above computations it is clear that the key ingredients in obtaining 
nonzero derivative F'([rf, r/“]) are sharp bounds for [a^] obtained from the cone 
condition and the transversality coming from dynamics, i.e. large in absolute 

value coefficient (W). 


6.4 Implementation notes. 

A short C++11 program that realizes part of a computer assisted proof of 
Theorem [T] is available from [Wj. All the inequalities, i.e. return times (43), 
shooting (441 and estimation on derivatives (46), regarding all six connecting 


orbits are checked within 0.7 second on a laptop type computer with Intel Core 
i7 2GHz processor. The program has been compiled and tested under Ubuntu 
OS with gcc-4.8.1 and gcc-4.9.1, OSX 10 with clang 6.0 and MS Windows 7 
with gcc-4.8.1. 

This program does not contain nonrigorous routines for finding good candi¬ 
dates for connecting orbits. They have been found by a simple bisection algo¬ 
rithm. The sizes of sets (see <m>) were adjusted by hand so that no subdivisions 
are necessary when computing derivatives F' n . 
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